Sorting of Chiral Microswimmers 
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Microscopic swimmers, e.g., chemotactic bacteria and cells, are capable of directed motion by ex- 
erting a force on their environment. For anisotropic microswimmers, e.g., bacteria, spermatozoa and 
many artificial active colloidal particles, a torque is also present leading to a circular motion with 
a well-defined chirality. Here, we demonstrate with numerical simulations how the chirality of this 
motion couples to chiral features present in the microswimmer environment. Levogyre and dextro- 
gyre microswimmers can be separated and selectively trapped in chtral flowers of ellipses. Patterned 
microchannels can be used as funnels to rectify the microswimmer motion, as sorters to separate 
microswimmers based on their linear and angular velocities, and as sieves to trap microswimmers 
with specific parameters. These phenomena can be scaled down to smaller microparticles as long as 
the Peclet number is maintained constant. 



Differently from simple Brownian particles, whose mo- 
tion is dominated by random thermal fluctuations, mi- 
croswimmers are capable of directed motion. Many 
microorganisms can actively explore their environment, 
e.g., E. coli bacteria and white blood cells Re- 
cently, several artificial microswimmers have also been 
demonstrated, e.g., self-diffusiophorctic Janus micropar- 
ticles ill. These microswimmers hold the promise of 
major applications in several fields, e.g., to deliver drugs 
within tissues, to localize pollutants in soils and to per- 
form tasks in lab-on-a-chip devices l3- 17j|. 

In order to perform active motion, a microswimmer 
must exert a dipole force on its environment 18, 3]. If 



this driving force acts along the line of motion, the mi- 
croswimmer will move along a straight line just perturbed 
by random Brownian fluctuations 6,]. However, often 
the microswimmer is anisotropic so that the driving force 
and the propulsion direction are not aligned [2^. This 
results in a torque acting on the microswimmer and in 
chiral active Brownian motion, i.e. motion along circular 
trajectories perturbed by Brownian fluctuations. There 
are many natural examples of chiral microswimmers. For 
example, E. coli bacteria (Fig. [T^) and spermatozoa un- 
dergo chiral active Brownian motion when moving near 
boundaries [21.-26]. Also, artificial microswimmers show 
characteristic chiral trajectories when they are asymmet- 
rical either by engineering or by chance. Moving down to 
the nanoscale, we can also consider molecular-sized chi- 
ral microswimmers, which can be obtained by joining a 
chiral molecule with a chiral propeller, e.g., a flagellum 
(Fig.Hb). 

Sorting microswimmers based on their swimming prop- 
erties, e.g. velocity, angular velocity and chirality, is 
of utmost importance for various branches of science 
and engineering. Genetically engineered bacteria can be 
sorted based on phenotypic variations of their motion [l[ . 
Velocity-based spermatozoa selection can be employed to 
enhance the success probability in the artificial fertiliza- 
tion techniques 113] • Considering the intrinsic variability 
of microfabrication techniques, the efficiency of artificial 



microswimmers for a specific task, e.g., drug-delivery or 
bioremediation, can be increased by selecting only the 
ones with the most appropriate swimming properties. 
Finally, the separation of levogyre and dextrogyre chi- 
ral molecules can be more effectively accomplished by 
chemically coupling them to a chiral propeller, sorting 
the resulting chiral microswimmer and finally detaching 
the propeller. This is important because, e.g., often only 
one specific chirality is needed by the chemical and phar- 
maceutical industry [281 ] and can hardly be achieved by 
mechanical means due to the extremely small Reynolds 
numbers [iol ]. 

In this article, we numerically demonstrate that chi- 
ral microswimmers can be sorted on the basis of their 
swimming properties by employing some simple static 
patterns in their environment. We show that a chiral 
flower, i.e. a chiral structure formed by some tilted el- 
lipses arranged in a circle, can trap microswimmers with 
a specific chirality and can, therefore, be used to sep- 
arate a racemic mixture. We also demonstrate that a 
patterned microchannel can be used as a funnel to rec- 
tify the motion of chiral microswimmers, as a sorter of 
microswimmers based on their linear and angular veloci- 
ties, and as a sieve to trap microswimmers with specific 
parameters. All these phenomena can be scaled down 
to smaller microparticles as long as the Peclet number 
is maintained constant. We studied microswimmers in 
two dimensions, but our results can be straightforwardly 
extended to the three-dimensional case. 



MODEL 

The motion of chiral microswimmers arises from 
the combined actions of random diffusion, an internal 
self-propelling force and a torque [l^. The position 
[x{t), y{t)] of a spherical particle with radius R undergoes 
Brownian diffusion with translational diffusion coefficient 



Dt = 



GnriR 
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TABLE I. Microswimmer parameters used in the simulations. From the radius R the rotational diffusion coefficient Dr and 
the translational diffusion coefficient Dt are obtained using Eqs. ^ and ^ respectively. The linear velocity v and the angular 
velocity uj are rescaled in order to maintain the Peclet number Pe = constant. 



R (pim) 


Dr (rad^s-') 


Dt (/ini^s-^) 


V (/ims 


n(rads"^) 


Pe 


1.00 


0.16 


0.22 


31.25 


±3.1 


142 


0.50 


1.32 


0.44 


125 


±12.6 


142 


0.25 


10.54 


0.88 


500 


±50.3 


142 



where fc^ is the Boltzmann constant, T is the temper- 
ature and 7] is the fluid viscosity. The particle self- 
propulsion results in a directed component of the motion, 
whose speed v we will assume to be constant and whose 
direction depends on the particle orientation (f{t). (p{t) 
undergoes rotational diffusion with rotational diffusion 
coefficient 



(2) 



In chiral microswimmers, (p{t) also rotates with angular 
velocity as a consequence of the torque acting on the 
particle. The sign of determines the microswimmer chi- 
rality. The resulting set of Langevin equations describing 
this motion is: 



V sin{ip{t)) dt + sfLD^ dW^ 



dip{t) 
dx{t) 

dy{t) = V cos((f{t)) dt + dWy 



(3) 



where W^p, Wx and Wy are independent Wiener pro- 
cesses. Equations ^ are then solved using some stan- 
dard finite-difference numerical methods [SO, 311. We 



have assumed that the microswimmers are spherical, e.g., 
active Janus microparticles @, but these equations 
can be straightforwardly generalized to elongated mi- 
croswimmers, e.g., nanorods and bacteria, by eniploying 
diffusion matrices instead of diffusion constants [32|. 

In patterned environments, the microswimmers will of- 
ten encounter obstacles 12[. In these interactions, when 
a microswimmer gets in contact with an obstacle, it 
landslides along the obstacle until its motion orientation 
points away from the obstacle. Numerically, this is im- 
plemented by suppressing the motion component perpen- 
dicular to the obstacle perimeter when pointing onwards 
the obstacle. 



HOMOGENEOUS ENVIROINMENTS 

The red line in Fig. [TJ; is a sample trajectory for a lev- 
ogyre microswimmer with R = 1.00 /xm (see Tab. |T] for 
the full list of parameters) moving in an environment 
free from obstacles. It bends counterclockwise track- 
ing almost circular trajectories just disturbed by Brow- 
nian fluctuations. Changing the chirality sign to dex- 
trogyre (black line in Fig. [TJ;), the trajectory behaves 



similarly but bending clockwise. This characteristic be- 
havior becomes clearer considering the ensemble average 
of many trajectories. In Fig.[Tf, we plot the average of 10^ 
trajectories starting at [a;(0), j/(0)] — [0/im, 0/im] with 
ip{0) = with levogyre (red line) and dextrogyre (black 
line) chirality; the average trajectories describe a spira 
mirabilis whose orientation depends on the motion chi- 
rality |2(| . The dimensions of the spire mirabilis define 
a length scale for the rotation of the swimmers, which is 
relevant in the presence of patterns in the environment. 
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FIG. 1. Chiral microswimmers. (a) Eschenchta coli bacteria 
perform a characteristic chiral motion in the proximity of a 
surface, (fe) Active chiral molecules can be obtained by chem- 
ically attaching a chiral molecule with a chiral propeller, e.g., 
a flagellum. (c — e) The trajectories of chiral levogyre (red) 
and dextrogyre (black) microswimmers with different radii 
[R = 1.00, 0.50 and 0.25 /im for (c), (d) and (e) respectively, 
see Tab. U for the other parameters) are qualitatively similar 
as long as the Peclet number is kept constant and the time 
is scaled accordingly {t = 10, 2.5 and 0.625 s for (c), (d) and 
(e) respectively). See also the supplementary movies 1, 2 and 
3 corresponding to (c), (d) and (e) respectively. (/) Aver- 
age of 10^ trajectories starting at [a;(0), y(0)] = [0 /im, /im] 
with v'(O) = for levogyre (red) and dextrogyre (black) mi- 
croswimmers as in (c). 



Similar results are obtained for smaller particles as long 
as the Peclet number is kept constant. This is achieved 
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by rescaling t, v and according to R~^. For example, in 
Figs. [TJi and e, the trajectories for R = 0.50 and 0.25 
respectively (see Tab. |T] for the fuU list of parameters) are 
qualitatively very similar to the ones for R = 1.00 /im 
(Fig.ttt). 

PATTERNED ENVIROINMENTS 
Chirality Separation 

The microswimmers can be selected on the basis of 
the sign of their motion chirality in the presence of some 
chiral patterns in the environment, e.g., an arrangement 
of tilted ellipses along a circle forming what we call a 
chiral flower. The shaded areas in Fig. [2^ depict a chi- 
ral flower where the ellipses are arranged along a circle 
with radius 11 /im. The dimensions of this structure are 
chosen to be comparable with the characteristic length 
scale defined by the spire mirabilis ([T]). The levogyre 
microswimmer (red line in Fig. [2^) is able to enter and 
exit the chiral flower without difficulty; however, a dex- 
trogyre microswimmer (black line in Fig. [2^) is trapped 
within the chiral flower. Fig.[2j) shows the radial position 
distribution of 10^ microswimmers after 10 s from when 
they are released from the center of the chiral flower: the 
levogyre microswimmers (red histogram) are outside the 
chiral flower, while the dextrogyre ones (black histogram) 
are trapped inside. 



(a) (b) 
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FIG. 2. Chiral microswimmers in a chiral environment, (a) 
10 s trajectories of a levogyre (red line) and dextrogyre (black 
line) microswimmer with R = 1.00 /xm (Tab. H]) in a chiral 
flower of ellipses (shaded areas). See also the supplementary 
movie 4. (b) Radial position distribution after 10 s of 1000 
levogyre (red histogram) and dextrogyre (black histogram) 
microswimmers starting from the center of the chiral flower 
at t = s. 

We can use the previous observation to devise a simple 
device to separate and trap microswimmers with differ- 
ent chiralities. As shown in Fig. [3l we use two chiral 
flowers with opposite chiralities enclosed in a box where 
the microswimmers can move freely. We start at time 
t = Os (Fig. [3^) with a racemic mixture placed inside 



each of the chiral flowers. Already at t = 10 s (Fig. [3]d), 
most of the levogyre (dextrogyre) microswimmers have 
escaped the right (left) chiral flower, while the oppo- 
site chirality ones remain trapped. In the subsequent 
time, the free microswimmers explore the space of the 
box util they get trapped by the corresponding chiral 
flower (Fig. [3j;-d). Almost all microswimmers are sta- 
bly trapped hy t = 500 s (Fig. [3^-f). We have verified 
that a similar behavior is obtained also for smaller mi- 
croswimmers {R = 0.50 and 0.25 ^i^Jl, Tab. H, ahhough 
for different parameters of the chiral flowers. 
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FIG. 3. Separation of levogyre and dextrogyre microswim- 
mers. (a) At t = 10 s, the microswimmers (7? = 1.00 /xm. 
Tab. lU are released inside two chiral flowers with opposite 
chirality. (6 — /) As time progresses, the levogyre (black sym- 
bols) microswimmers are trapped in the right chiral flower and 
the dextrogyre (red symbols) ones in the left chiral flower. See 
also the supplementary movie 5. 



Sorting by Velocity 

A patterned microchannel can be used to sort the mi- 
croswimmers on the basis of their linear velocity. We 
show an example of such microchannel in Fig. 2^ (grey 
areas). It is formed by a series of elliptical structures 
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inserted on the channel wahs at an angle and slightly 
shifted between the top and bottom sides. In this way 
the channel is itself chiral. Chiral microswimmers placed 
at position x = /im at i = s behave differently depend- 
ing on their linear velocity. For example, a.t t = 1000 s, 
levogyre microswimmers with v = 70iim./s (red sym- 
bols in Fig. |4^) have propagated towards the right sev- 
eral hundreds microns, while slower microswimmers with 
V — AO fim/s (black symbols in Fig. |4^) are trapped near 
the initial position; both microswimmers have the same 
chirality and angular velocity f2 = -1-3.1 rad/s. The pat- 
terned microchannel works as a funnel to rectify the mo- 
tion of chiral microswimmers, as demonstrated from the 
fact that the microswimmers move towards the right. 
The microswimmer motion in the channel can be thought 
of as a ratchet driven by the microswimmers own self- 
propulsion ||33j]. The microchannel works also as a sieve 
to trap microswimmers with specific parameters, e.g., in 
Fig.|4]only microswimmers with v > 40/im/s can propa- 
gate to the right. Changing the structure of the channel 
it is possible to change the velocity threshold below which 
the microswimmers start propagating. Finally, it is also 
possible to use the channel as a microswimmer sorter 
based on their linear, as shown in Figs. IDd-c. When we 
place at t = s a mixture of particles with t; = 40, 50, 60 
and 70 /xm /s at position x = /xm / s (Fig. [4)d) , we observe 
that the probability distributions of particles with differ- 
ent speed separate over time (Figs.|4l:-e). Faster particles 
propagate further along the channel while the slower par- 
ticles are held back. The separation efficiency increases 
with time as the particles propagate further along the 
channel. 



Sorting by Angular Velocity 

The patterned microchannel introduced in the previ- 
ous subsection can also be used to sort particles accord- 
ing to their angular velocity. In Fig. we show that 
microswimmers with v = 40//m/s and Q = 3.1 rad/s 
(black symbols) are trapped at their initial position af- 
ter 10000s, while microswimmers with ft — 2.2 rad/s 
(red symbols) can propagate. Over time, it is possible 
to separate particles with smaller differences in angu- 
lar velocity, as shown by the histograms in Fig. |4l3-e for 
n = 2.2, 2.5, 2.8 and 3.1 /im/s. 



CONCLUSIONS AND OUTLOOK 
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FIG. 4. Linear velocity based sorting in a microchannel. (a) 
Patterned microchannel (grey areas) and position of levogyre 
1.00 /im microswimmers with ?; = 40 (black symbols) and 
70 /xm/s (red symbols) 1000 s after they have been released 
from position x = 0/im (other parameters as in Tab.|l]). See 
also the supplementary movie 6. (6 — e) Histograms at vari- 
ous times of 1000 particles with v — 40, 50, 60 and 70 /xm/s 
released at t = Os from x — 0/im. 



ods, which only rely on static environmental structures 
and on the microswimmers' own self-propulsion, with- 
out requiring moving parts, fluid flows or external forces. 
These techniques can be exploited to separate naturally 
occurring chiral active particles, e.g., bacteria and sper- 
matozoa, and also to separate chiral passive particles, 
e.g., chiral molecules, by coupling them to some chiral 
motors. 
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We have studied the sorting of chiral microswimmers 
on the basis of the chirality, linear velocity and angu- 
lar velocity of their motion. This is obtained by using 
an environment with chiral patterns, which can be read- 
ily fabricated by standard microfabrication techniques. 
The main advantage lays in the simplicity of these meth- 
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